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Abstract 

Background: Panel studies have shown adverse effects of air pollution from biomass burning on children's health. This 
Study estimated the effect of current levels of outdoor air pollution in the Amazonian dry season on peak expiratory flow 
(PEF). 

Methods: A panel study with 234 schoolchildren from 6 to 15 years old living in the municipality of Tangara da Serra, Brazil 
was conducted. PEF was measured daily in the dry season in 2008. IVlixed-effects models and unified modelling repeated for 
every child were applied. Time trends, temperature, humidity, and subject characteristics were regarded. Inhalable 
particulate matter (PMio), fine particulate matter (PIVI2.5), and black carbon (BC) effects were evaluated based on 24-hour 
exposure lagged by 1 to 5 days and the averages of 2 or 3 days. Polynomial distributed lag models (PDLIVl) were also 
applied. 

Results:The analyses revealed reductions in PEF for PMio and PMj.s increases of 1 0 |jg/m^ and 1 |ig/m^ for BC. For PMio, the 
reductions varied from 0.15 (confidence interval (Cl)95%: -0.29; -0.01) to 0.25 l/min (CI95%: -0.40; -0.10). For PIVI2.5, they 
ranged from 0.46 {CI95%: -0.86 to -0.06) to 0.54 l/min (CI95%:-0.95; -0.14). As for BC, the reduction was approximately 
1.40 l/min. In relation to PDLIVl, adverse effects were noticed in models based on the exposure on the current day through 
the previous 3 days (PDLM 0-3) and on the current day through the previous 5 days (PDLIVl 0-5), specially for PMio- For all 
children, for PDLIVl 0-5 the global effect was important for PMio, with PEF reduction of 0.31 l/min (CI95%: -0.56; -0.05). 
Also, reductions in lags 3 and 4 were observed. These associations were stronger for children between 6 and 8 years old. 

Conclusion: Reductions in PEF were associated with air pollution, mainly for lagged exposures of 3 to 5 days and for 
younger children. 
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Introduction 

Several studies have shown adverse health effects from air 
pollution in many parts of the world [l]-[3]. Panel studies with 
children and teenagers have revealed important associations 
between air pollution and episodes of respiratory symptoms or 
lung function [4]-[8]. Such effects were found in both asthmatic 
and healthy subjects aged between 6 and 18 years old. Because 
children are stiU developing physiologically, exposure to air 
pollutants is a major risk factor for health and may have 



consequences in adulthood [9]. However, air pollution studies in 
healthy children are scarce. 

In Brazil, the health effects of air pollution from biomass 
burning have not been thoroughly studied, especially in the 
Amazon region. The Brazilian Amazon dry season is the most 
critical period for biomass burning. From June to October, levels 
of particulate matter are usually above the World Health 
Organisation guidelines [2]. 

Some studies undertaken in the Amazon region have shown 
adverse health effects from particulate matter due to biomass 
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Figure 1. Geographic location of the municipality of Tangara da Serra, Mato Grosso State. 

doi:! 0.1 371/journal.pone.01 041 77.g001 



burning, for instance, increases in emergency room visits, 
outpatient visits, and liospital admissions as well as decreases in 
peak expiratory flow (PEF) and increases in the fi'cquency of 
micronuclei in oral epithelial cells [10]-[15]. These effects were 
found in children and the elderly in municipalities around the Arc 
of Deforestation region. 

The Brazilian Amazon covers an area of approximately 
5,000,000 km and 9 Brazilian states. Over 25 million inhabitants 
are changing the land use [16]. The majority of forest fire hotspots 
takes place in the Arc of Deforestation [17], whilst the state of 
Mato Grosso usually shows the worst records. Ignotti et al. [17] 
compared the indicators of morbidity and mortality caused by 
respiratory diseases in municipalities of Mato Grosso. Alta Floresta 
and Tangara da Serra were identified as priority areas for assessing 
the air pollution health effects of biomass burning in the Brazilian 
Amazon. 

A panel study of schoolchildren carried out in the municipality 
of Alta Floresta found adverse effects of particulate matter with an 
aerodynamic diameter less than 2.5 |J,m (PM2.5) on PEF [10]. The 
results were stratified by the time of day the children attended 
school, morning or afternoon. Stronger associations were detected 
in the afternoon group. However, only single-lag exposures were 
estimated, namely, levels of PM2.5 on the current day or lagged by 
1 or 2 days and the averages of 0- to 1-day lags, 1- to 2-day lags, 
and 0- to 2-day lags. 



A similar panel study was carried out in Tangara da Serra in 
2008. Its results are reported in this article. The primary aim of 
this study was to estimate the effect of current air pollution levels in 
the Amazonian dry season on schoolchUdren's PEF. Three 
separate secondary aims were investigated, namely: 

(i) To compare two methodological approaches in the analysis 
of repeated measures data: mixed-effects models (MEM) and 
a unified model to be applied for every child; 

(ii) To assess the role played by age in explaining the difference 
between the morning and afternoon groups of children that 
Jacobson et al. [10] uncovered. In Brazil, most primary 
schools adopt the policy of two periods of study classes, 
morning and afternoon. Further, younger children attend 
school in the afternoon and older children go to school in 
the morning. Therefore, age could explain the difference 
observed in the previous study; 

(iii) To extend the lag structures that Jacobson et al. [10] 
explored but to also use a polynomial distributed-lag model 
(PDLM). In the Amazon, particles emitted by biomass 
burning remain in the atmosphere for one wcc-k or so [18]. 
Therefore, to evaluate the cumulative effects based on 
polynomial distributed-lag models seems to be more 
appropriate. Moreover, few panel studies have used PDLM 
to estimate air pollution effects, especially those related to 
children's PEF [19]-[21]. 
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Material and Methods 

Study Design and Variables 

The study was conducted in the Tangara da Serra municipahty, 
in the southwest of Mato Grosso State, in a transition area of the 
Amazon and Cerrado biomes (Figure I). Its total population was 
83,431 inhabitants (population density was 7.32 inhab./km% and 
the population from 6 to 15 years old was 14,398 inhabitants 
(17.3%) [16]. Its main economic activities are agriculture and 
livestock. Despite its increasing production of sugarcane, the 
harvest is still manual and occ:urs during the dry season [22] . The 
pre-harvest burning of sugarcane and the plumes of biomass 
burnings from neighbouring regions are the main sources of air 
pollution in the city. 

In 2008, a panel study of schoolchildren from a local public 
school was carried out (Lat 14°37'59.05" S; Long 57°3r25.34" 
W). The school had 875 students from 6 to 15 years old attending 
school in the morning or in the afternoon. The study sample 
comprised 234 schoolchildren who were followed up daily, except 
weekends and holidays, from August to November. Children were 
randomly selected according by age group (between 6 and 15) and 
school schedule. 

Ethical appr()\'al was obtained from the UNEMAT (Mato 
Grosso State University) Ethics Committee, and a consent form to 
participate in the study was signed by the children's parents or 
legal guardians. The protocol described by the American Thoracic 
Society [23] for PEF exams was adopted throughout the study. 
After careful training, nursing undergraduate students collected 
PEF data under the supervision of a senior pulmonologist and a 
senior nurse. 

Three daily PEF measures were sequentially taken. To avoid 
disruption of the school routine, groups of 4 or 5 children 
performed the tests and returned to class afterwards. The device 
used was the Mini-Wright Peak Flow Meter (Clement Clarke 
International Ltd, United ICingdom). The study outcome was the 
maximum of the three sequential PEF measures. 

Weight and height were measured at the study's onset, using a 
Welmy anthropometric mechanical scale with 150 kg maximum 
capacity. Face-to-face questionnaires comprising health conditions 
and exposure to tobacco smoking at home were answered by the 
children's guardians/parents. To assess asthma symptoms, the 
questionnaire included phase I of the International Study of 
Asthma and Allergies in Childhood [24] . 

Exposure to air pollution was measured during the entire study 
period. The devices were placed at the UNEMAT campus in 
Tangara da Serra. The selected pupils lived within a 5 km radius 
of the campus and a maximum of 1 km radius of the study school. 
PMio concentrations (particulate matter with an aerodynamic 
diameter less than 10 Jtm) were measured continuously using the 
Thermo Electron Corporation (DATARAM) device. Moreover, 
the Aethalometer was used to measure blac:k carbon (BC) levels, 
and stacked-filter units (SFU) were also used to measure PM2.5. 
The following were regarded as exposures to PMio and BC: (i) the 
24-hour mean of the previous day (lag 1) until 5 days before (lag 5); 
and cumulative exposures, such as the average of the daily means 
of the same day and the previous one (lags 0-1), of the previous 
two days (lags 1-2), and of the same day and the two previous days 
(lags 0-2). For PM2.5, there were a few missing data; therefore, to 
minimise loss of information, cumulative exposures were not 
investigated. 

Measures of relative humidity and temperature were provided 
by INMET (National Institute of Meteorology, Brazil). Levels of 
relative humidity and mean temperature were regarded in all 
statistical models. 



Statistical Analysis 

The study database included variables that vary over time, for 
example, the outcome, the exposure, the weather, and outdoor 
leisure time, but also time-invariant covariates, such as age, sex, 
passive smoking, asthma, weight, height, and body mass index 
(BMI). Despite the data collection from August to November, this 
study presents the results for the dry season period only, from 
August to October. 

Five types of statistical analysis were implemented in order of 
complexity, namely: (i) separate data analysis for each child; (ii) a 
mixed-effects core model regarding climate and individual 
covariates; (iii) air pollution effect estimation for single lags; (iv) 
air pollution effect estimation for PDLM; (v) and a sensitivity 
analysis. 

For item (i), the same modelling steps were applied to the PEF 
time series for each child. Several strategies were examined for the 
adjustment of time trends, temperature and humidity, such as 
linear, quadratic or parametric splines with one knot or two knots, 
as well as linear or quadratic polynomials. The modelling 
proceeded including a 1*' order auto regressive term, later 
estimated for each child. To summarise the results, averages were 
calculated for the Akaike Information Criterion (AIC), the 
estimate of the P' order autoregressive term ), and the set of 
PMio lag 1 effects. Positive (p<0.05), negative (p<0.05), or null 
effects (p>0.05) of PMjo lag 1 were classified into three groups. 
The unified model that presented the smallest AIC average was 
selected. Finally, tables of counts (proportions) or averages by sex, 
age, passive smoking, presence of asthma, BMI, weight, and height 
allowed comparison of the three groups. 

Because this was a study of repeated measures, MEMs were also 
appropriate for the statistical analysis. Regarding this approach, 
the model adjustment considered a common autoregressive lag 1 
(AR(1)) correlation structure whereas the variance function 
comprised three time-invariant covariates, namely, BMI in four 
categories, age in three, and asthma status. Four adjustments 
based on polynomials of time were evaluated for the outcome 
long-term trend: a cubic parametric spline with one knot; a 
quadratic parametric spline with two equally spaced knots; a 
quadratic parametric spline with one knot; and a quadratic 
polynomial of time [25]. The above considered time to be centred 
in the middle of the study period as well as random coefficients. 
The same approach was applied for temperature and relative 
humidity, but also linear effects were explored. Furthermore, 
exposures on the same day, the previous day, and the day before 
that were examined. Finally, height, weight, BMI, age, sex, 
ISAAC asthma status, passive smoking, and outdoor leisure time, 
with fix(-d (-ffects, were regarded in the models. Because of the 
number of model comparisons, the AIC and the deviance test were 
used to select the most parsimonious model, which then was 
scrutinised with regard to residual diagnostics to confirm the 
normality assumptions of the random effects. Jacobson et al. [10] 
describe the details of this approach. 

Regarding the choice of core model based on the modelling 
steps mentioned above, single-lag and PDLM were applied to 
estimate the effects. Single-lag models were fitted for all exposures 
mentioned previously. PDLM estimates a fixed number of lagged 
effects constrained to a given choice of polynomial as well as the 
overall sum of effects [26]. In this study, quadratic polynomials of 
the effects lagged up to 3 (0-3) or 5 (0-5) days were calculated. 

Sensitivity analyses were performed to assess the robustness of 
the results. For instance, measures of PM2.5 were missing for 25% 
of the study days. Therefore, missing values of PM2.5 were 
replaced by the PM2.5 average on the previous day and the day 
after. The effect of this exposure on PEF was estimated using single 
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lags and PDLM. Because of the methods used in PDLM, we 
evaluated the effects of PM2.5 only considering levels after 
imputation. Furthermore, all of the analyses were repeated 
excluding outliers, the definition of which was standardised 
residuals greater than four in absolute value. We also explored 
the presence of a threshold effect of PMio using a linear 
parametric spline. Each exposure (lag 1, lag 2, lag 3, lag 4, lag 
5, and cumulative lag 1-5) was centred at different thresholds, 
namely, 20 Hg/m^, 30 \lg/m^, 40 \lg/m^, 50 Hg/m^, 60 \ig/m^, 
70 |J.g/m'\ and 80 |J.g/m\ 

All analyses were performed using R 2.15.1 [27] and the library 
nlme [28]. 

Results 

Of the 234 students selected for the sample, 7 left school during 
the study and were considered lost to follow-up. A group of 7 
students were excluded from the regression analysis because their 
guardians/parents did not answer the questionnaire. The total 
number of study days was 53. 

Table 1 displays the summary statistics of individual character- 
istics (sex, age, weight, height, asthma, and passive smoking), PEF 
measures, and weather and air pollution variables. Forty-seven 
percent of the 227 students were boys, the average age was 10.3 
years old, the average weight was 36.2 kg (17 to 82.2 kg), and the 
average height was 1.40 m (1.11 to 1.73 m). According to the 
ISAAC criteria, 18% of the students (n = 39) reported asthma 
symptoms (21 girls and 18 boys). Passive smoking at home was 
reported by 33% of the students' parents (n= 72). Overall mean 
PEF was 289.9 1/m and it ranged from 70 to 780 1/min. 

The average relative humidity and the daily mean temperature 
were 57% (29.4 to 90.5%) and 26°C (16.5 to 29.9°C), respectively. 
The averages of the daily mean PMio, PM2.5, and BC were 
62.7 ng/m^ (12.4 to 164.5 Hg/m'), 19.6 |J,g/m^ (4.4 to 55.9 ng/ 
m^), 0.998 ng/m'' (0.271 to 2.146 \ig/m% respectively (Table 1). 
Typically in this part of Brazil in October and November, 
humidity levels increase and air pollution levels decrease because 
of the beginning of the rainy season, whereas temperature levels 
are somewhat constant. Figure 2 shows the time series for 
temperature, humidity, PMio, PM2.5, and BC. 

In this section, the effects reported are in litres per minute and 
are for increases of 10 (tg/m^ for PMip and PM2.5, and 1 Jig/m^ 
for BC. 

Results for the Unified Model Approach for Each Child 

Table 2 shows the summary statistics of the unified model 
approach. Regarding AIC, l'"' order autocorrelation (Jp ), PMio lag 
1 effect averages, and counts of individuals with positive or 
negative effects of PMio lag 1 are displayed. The most 
parsimonious 1*' order autoregressive model consisted of a 
quadratic parametric spline with two equally spaced knots for 
the long-term trend, a quadratic polynomial for humidity lagged 
by two days, and temperature lagged by two days. The 0 average 
was 0.0232 (sd 0.246; -0.812 to 0.549), and for 34 children the 
fitted models presented ^ larger than zero and p<0.05. With 
respect to PMio l^ig 1, 12 students presented positive effects (p< 
0.05), 13 students presented adverse effects (p<0.05), and 202 
presented no effects (p>0.05). 

Table 3 summarizes the results of the AR(1) estimated model 
according to individuals characteristics. Averages of weight, 
height, BMI, and age were calculated according to these groups 
(positive, negative, and non-significant effect) and compared using 
the Kruskal-Wallis test. Of the four individual covariates, only age 
showed differences (p<0.05). Similarly, proportions of sex. 



asthma, and passive smoking were compared across groups using 
Fisher's exact test. There were no proportion differences. Of the 
individual covariates investigated, only age affected the magnitude 
of PMio effects (Table 3). 

Moreover, we computed combined-effect estimates of the 
exposure to PMio levels lagged by between 1 and 5 days, using 
averages weighed by the inverse of their variances. The most 
adverse effect found was for lag 4. For PMio lag 4, there was a 
reduction in PEF equal to 0.16 (95% confidence interval(CI): 
—0.31; —0.02). The other results are not shown. 

Results for the Mixed-Effects Models 

MEMs were applied in the next step of the data analysis. The 
choice of core model consisted of long-term trend adjusted by a 
quadratic parametric spline with two equally spaced knots and 
random coefficients, a quadratic polynomial cur\'e for humidity 
lagged by two days with random coefficients, temperature lagged 
by two days and a random coefficient, and fixed-effects covariates 
such as BMI, age, sex, asthma status, and exposure to air pollution; 
it also regarded an AR(1) and sixteen different subject-error 
variances depending on the subject status with respect to asthma, 
age, and BMI. Table 4 summarises the adjusted mixed models. 

Single-Lag Effects 

Table 5 displays the results of the models for all students and 
stratified by three age groups and considering single-lag exposures. 
This approach revealed some adverse effects of air pollution on 
PEF. Taking into account all children: (i) for PMio, the results 
showed reductions of 0.25 (95%CI: -0.40; -0.10), 0.20 (95%CI: 
-0.32; -0.07), and 0.15 (95%CI: -0.29; -0.01) in lags 3, 4, and 
5, respectively; (ii) for PM2.5 lag 4, there was a reduction of 0.54 
(95%CI: -0.95; -0.14). However, after PM2.5 imputation, the 
models revealed associations not only with lag 4 (—0.50; 95%CI: 
-0.89; -0.12) but also with lag 3 (-0.46; 95%CI: -0.86; -0.06); 
(iii) for BC, the reduction was roughly 1.40 (95%CI: -2.50; 
-0.29) in both lags 4 and 5. 

The results were stratified by age: 6 to 8 years old, 9 to 11, and 
12 to 15. Children aged 6-8 were the most susceptible group. The 
results regarding this group were: (i) for PMio, the reduction varied 
from 0.28 (95%CI: -0.50; -0.06) to 0.43 (95%CI: -0.65; -0.20) 
in lags 2, 3, 4, and 5, and there was also a reduction of 0.30 
(95%CI: —0.54; —0.06) for the cumulative exposure in lag 1-2; (ii) 
for PM2.5, important reductions were observed for lag 2 (—0.91; 
95%CI: -1.65; -0.16) and lag 4 (-0.93; 95%CI: -1.55; -0.32). 
After imputation, the adverse effect of PM2.5 was noticed only for 
lag 4 (-0.87; 95%CI: - 1.47; -0.28); (iii) for BC, lagged exposures 
of 4 and 5 days reduced PEF in this age group. Further, for 
children aged 9-1 1, the adverse effect was observed only for BC. 
For lag 5, there was a reduction of 1.88 (95%CI: -3.64; -0.12). 
No associations were found for children aged 12-15. 

Polynomial Distributed-lag Effects 

Table 6 presents the results of the PDLM 0-3 for all children 
and stratified by age. Some adverse effects were revealed under 
this approach. For PMio, the overall effect was not significant for 
an children; however, lag 3 showed a reduction of 0.25 (95%CI: 
— 0.43; —0.08). There was also an overall adverse effect of PMio 
for children 6-8 (-0.34; 95%CI: -0.66; -0.01), particularly for 
lag 3 (-0.35; 95%CI: -0.62; -0.08). For PM2.5 lag 3, the group 
of all children suffered a reduction of 0.41 (95%CI: -0.81; - 
0.01). There were no adverse effects for BC. Likewise, no evidence 
of adverse effects was found for children aged 9-15. 

The PDLM 0-5 effect estimates revealed stronger evidence of 
adverse effects (Table 7). Regarding PMio for all children, the 
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Figure 2. Daily temperature, Kiumidity and air pollutant readings. Tangara da Serra, MT, Brazil-2008. 
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Table 3. Characteristics of the schoolchildren according to the results of the AR(1) estimated model. 





Individual CKiaracteristlcs 


Change in PEF - PM^o ^-ag^ 






Positive and Significant (n = 12) 


Negative and Significant (n = 13) 


Non-Significant (n = 202) 


Sex (%) 


Girls 


66.7 (n = 8) 


46.2 (n = 6) 


53.0 (n = 107) 


Boys 


33.3 (n = 4) 


53.8 (n = 7) 


47.0 (n = 95) 


Asthma Status (%) 


Yes 


8.3 (n = 1 ) 


23.1 (n = 3) 


18.1 (n = 35) 


No 


91.7 (n = 11) 


76.9 (n = lC) 


81.9 (n = 158) 


Passive Smoke {%) 


Yes 


33.3 (n=4) 


30.8 (n = 4) 


33.2 (n = 64) 


No 


66.7 (n = 8) 


69.2 (n = 9) 


66.8 (n = 129) 


Age 


Average (SD) 


12 (2.17) 


9.62 (2.43) 


10.23 (2.55) 


Weight 


Average (SD) 


43.57 (15.01) 


35.26 (11.42) 


35.92 (13.40) 


Height 


Average (SD) 


149.3 (14.16) 


138.69 (11.93) 


139.71 (15.28) 


BMi 


Average (SD) 


19.05 (3.87) 


1 7.97 (4.03) 


17.79 (3.50) 



Tangara da Serra-MT, Brazil, 2008. 
doi:1 0.1 371 /journal.pone.Ol 041 77.t003 



overaU effect was -0.31 (95%CI: -0.56; -0.05). Additionally, 
reductions of 0.09 (95%CI: -0.16; -0.01) and 0.10 (95%CI: 
— 0.15; —0.05) in lags 3 and 4, respectively, were estimated. For 
children aged 6-8, there were strong effects for lags 2, 3, and 4. 
The overall effect for this age group indicated a reduction of 0.52 
(95%CI: -0.92; -0.12) caused by PMk,. 

Moreover, for all children, the PDLM 0-5 effect estimates of 
PM2.5 showed reductions for lag 2 (-0.21; 95%CI: -0.42; 
-0.003), lag 3 (-0.27; 95%CI: -0.48; -0.06), and lag 4 (-0.19; 
95%CI: -0.36; -0.02), but there was no overaU effect. 
Additionally, for children aged 6-8, the results revealed reductions 
for lags 3 and 4. 



BC presented associations with PEF. For all children, PDLM 0— 
5 revealed a reduction only for lag 4 (-0.73; 95%CI: -1.22; 
— 0.24). Furthermore, for children aged 6-8, there were reductions 
of 0.97 (95%CI: -1.74; -0.20) and 1.77 (95%CI: -3.26; -0.27) 
in lags 4 and 5, respectively. 

Sensitivity Analysis 

Exclusion of outliers hardly affected the results. However, other 
exposures showed adverse effects, for instance: (i) For single-lag 
models, PMio lag 3 presented a negative association with PEF for 
children aged 6-8, and PM2.5 lagged by 3 days for all children; (ii) 
for PDLM, a reduction was observed with PM2.5 lag 2 in children 
aged 6-8. 



Table 4. Summary adjustment strategies of the mixed-effect models. 



Model Adjustment* (n = 220) ^IC PHi(AR(1)) Change in PEF ••(iC9S%) 



Time trend 


Humidity (lag 2) 


Temperature (lag 2) 










Linear 


Linear 


Linear 


90301.19 


0.332 


-0.020 (- 


-0.165; 0.125) 


Linear 


Quadratic Polynomial 


Linear 


90292.05 


0.330 


-0.021 (- 


-0.165; 0.123) 


Linear Spline^ 


Linear 


Linear 


90207.98 


0.291 


-0.016 (- 


-0.157; 0.125) 


Quadratic Polynomial 


Quadratic Polynomial 


Linear 


90217.62 


0.295 


-0.026 (- 


-0.168; 0.115) 


Quadratic Spline^ (1 Knot) 


Quadratic Polynomial 


Linear 


90147.36 


0.261 


-0.048 (- 


-0.188; 0.092) 


Quadratic Spline^ (2 Knot) 


Quadratic Polyriomial 


Linear 


90104.27 


0.238 


-O108 (- 


-0.251; 0.035) 


Cubic Spline^ (1 Knot) 


Quadratic Polynomial 


Linear 


90120.02 


0.242 


-0.114 (- 


-0.259; 0.032) 



Tangara da Serra-MT, Brazil, 2008. 

*Model adjusted for age, BMI, gender, and asthma status; random coefficients for the intercept, time trend, humidity, and temperature; variance function of the random 

error included age, BMI, and asthma status. 

** Change in PEF for an Increase of 10 fig/m^ In PMiolagl. 

^Parametric Splines. 

doi:l 0.1 371/journal.pone.Ol 041 77.t004 
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Figure 3. Estimated changes in peak expiratory flow (in l/min) for an increase of 10 {ig/m^ in PM^oi above the thresholds for all 
children. Tangara da Serra-MT, Brazil - 2008. 
doi:1 0.1 371/journal.pone.01 041 77.g003 



We found evidence for PMio thresholds. Under the linear 
parametric spline approach, there were two slopes, one for the 
effect estimate below the threshold and one for the effect above it. 
Concerning the former, there were no important ffndings, except 
for lags of PMio such as: (i) lag 2: there was a reduction of 0.70 
(95%CI: -1.30; -0.99) for exposures below 50 ^ig/m'^; (ii) lag 3: 
the effects of exposures below 60 |J,g/m', 70 \ig/ m\ and 80 \ig/ 
m^ were reductions of 0.49 (95%CI: -0.93; -0.05), 0.45 (95%CI: 
-0.83; -0.07), and 0.40 (95%CI: -0.72; -0.08), respectively; (iii) 
lag 1-5: reductions for the exposures below 60 |ig/m' (-0.63; 
95%CI: -1.18; -0.08), 70 |ig/m'' (-0.55; 95%CI: -0.99; 
-0.11), and 80 [ig/rn (-0.46; 95%CI: -0.84; -0.09). 

Figure 3 presents a bubble plot for the estimated effects above 
the thresholds for all exposures. The size of the bubbles is 
proportional to the number of days above the thresholds. Negative 
associations between PEF and PMio were noticed for lags 3, 4, 5, 
and lag 1-5. More specifically, for lag 3 the associations were 
important for the following thresholds: 20 Hg/m"" (-0.25; 95%CI: 



-0.40; -0.10); 30 ng/m' (-0.26; 95%CI: -0.41; -0.11); and 
40 |tg/m'^ (-0.23; 95%CI: -0.40; -0.06). For lag 4, the effects 
were considerable for the same thresholds, varying from 
-0.20 ng/m'^ to -0.21 |ig/m'^ (95%CI: -0.35, -0.08), as well 
as for 50 Hg/m'^ (-0.18; 95%CI: -0.34; -0.02). Concerning lag 
5, there were reductions for all thresholds: (i) 20 |tg/m'^, —0.15 
(95%CI: -0.30; -0.01); (ii) 30 |Xg/m^ -0.20 (95%CI: -0.35; 
-0.04); (iii) 40 \ig/m\ -0.19 (95%CI: -0.36; -0.03); (iv) 50 [ig/ 
m\ -0.23 (95%CI: -0.42; -0.05); (v) 60 ng/m\ -0.26 (95%CI: 
-0.47; -0.05); (vi) 70 |J.g/m'\ -0.28 (95%CI: -0.53; -0.03); (vii) 
80 \ig/m\ -0.32 (95%CI: -0.62; -0.03). For the cumulative 
exposure (lagl-5), the reductions were associated with the 
following thresholds: 20 \lg/rn (-0.37; 95%CI: -0.57; -0.16); 
30 ng/m'^ (-0.35; 95%CI: -0.55; -0.14); and 40 |.ig/m'^ (-0.37; 
95%CI: -0.61; -0.13). 
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Discussion 

This study's results suggest that air pollution from biomass 
burning may be a respiratory health risk factor for schoolchildren 
aged 6 to 15 years old living in the Brazilian Amazon. This 
association is consistent with the previous panel study undertaken 
in the municipality of Alta Floresta [10]. The current study 
extended previous analyses by including PDLM as well as adding 
the effects of BC and PMio exposures. Further, we explored two 
methods for analysing panel studies: MEM and univariate time 
series modelling, which was applied to every child. Although the 
combined-eflFect estimates obtained using the univariate approach 
and those from MEM were similar, their precision was different. 
However, it was particularly important to show that some students 
were more vulnerable to air pollution than were others. 

Our findings corroborate the results of panel studies describing 
a linear effect of air pollution on PEE [5], [8], [29]. No association 
was found with PM2.5 or BC in a rural panel study undertaken in 
the southeast of Brazil to investigate the effects of pre-harvest 
sugarcane burning on PEE [30]. Systematic reviews of panel 
studies demonstrated negative pooled effects of PMio and PM2.5 
on PEE [31], [32], [33]. 

Although we found associations for all children, the effects were 
stronger after stratification by age, with PEE decrements for the 
youngest group, confirming that younger children are more 
susceptible to air pollution effects [1]. In the Alta Eloresta study 
[10], a negative association of PM2.5 with PEE was found for 
children who studied in the afternoon. However, the authors noted 
that the majority of children in this group were from 6 to 9 years 
old. Our results confirm the hypothesis that age is most likely to 
explain this finding. 

A multitude of factors are related to children's PEF. In the 
present study, PEF measurements were associated with temper- 
ature, humidity, BMI, age, sex, and asthma status. Weather 
conditions can influence children's health [34] but also interfere 
with pollutants in the atmosphere [35]. Other environmental 
factors could affect PEF, for instance, time spent outdoors and 
passive smoking. However, in this study these variables were not 
important in the modelling, most likely because in the dry s(;as()n 
children spend most of their time outdoors when they are not at 
school because of lack of leisure indoors as well as the high 
temperatures. 

Most adverse effects found in this study were lagged by 3, 4, and 
5 days. Lagged effects as well as (■umulativ(- (■ff(;cts are expected in 
this region not only because of the characteristics of Amazon 
biomass burning (every day during the dry season) but also 
because of meteorological factors, which allow for pollutants to 
remain in the atmosphere for long periods of time. The PDLM 
approach allowed exploring the effect of air pollutants on PEF 
distributed over time and also the overall effect up to 3 or 5 lagged 
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